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ABSTRACT: The present paper is dealt with the study of the exact non-vacuum perfect fluid solutions of five 
dimensional Bianchi-type-I space-time in the metric version of /(/?) gravity considering stiff matter to obtain 
energy density and pressure of the universe on the lines of M. Sharif and M. Farasat Shamir (2010). In 
particular, we obtain two exact solutions which correspond to two models of the universe in higher five 
dimensional Bianchi type-I space-time. The function f(R) are also evaluated for both the models. Finally the 
physical properties of these models have been discussed. 
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I. INTRODUCTION 

S. N. Pandey (2008) has developed a higher order theory of gravitation based on a Lagrangian density 
consisting of a polynomial of scalar curvature R to obtain gravitational wave equations conformally flat. 
Recently S. N. Pandey and B. K. Sinha (2009) have studied spherically symmetric metric in the field equations 
of higher order theory of gravitation which is obtained by modifying Einstein's field equations in general 
relativity theory. Thus Einstein's general relativity is nothing but the particular case of this higher order 
modified theory i.e, f(R) theory of gravity. 

Many authors have obtained spherical, cylindrical and plane symmetric solutions in f (R) gravity 
[Hollenstein and Lobo (2008), Azadi et al.(2008), M. Sharif and M. Farasat Shamir (2009), etc.]. Sotiriou and 
Faraoni (2010) presented some important aspects of f(R) theory of gravity in Metric, Palatini and Metric- 
Affine formalisms. 

Recently, in the paper [1], M. Sharif and M. Farasat Shamir (2010) have studied the non-vacuum 
perfect fluid solutions of four dimensional Bianchi type-I space-time in the framework of metric f (R) gravity. 

For this purpose they have considered stiff matter to find energy density and pressure of the universe and found 
two exact solutions which correspond to two models of the universe and observed that the first solution yields a 
singular model while the second solution gives a non-singular model. Finally the physical behavior of these 
models has been discussed using some physical quantities. The function f ( R) of the Ricci scalar is evaluated 

also by them. It is noted that the work of M. Sharif and M. Farasat Shamir (2010) is in the framework of four 
dimensional Bianchi type-I space-time. In the present paper we wish to investigate the similar situation in higher 
five dimensional Bianchi type-I space-time. 

In a past few years there have been many attempts to construct a unified field theory based on the idea 
of multidimensional space-time. The idea that space-time should be extended from four to higher five dimension 
was introduced by Kaluza and Klein (1921, 26) to unify gravity and electromagnetism. Several aspects of five 
dimensional space-time have been studied in different theories by many authors [Wesson (1983, 84), Reddy 
D.RK. (1999), 

Khadekar et al. (2001), Ghosh and Dadhich (2001), Adhao (1994), Thengane (2000), Ambatkar (20002), 
Jumale(2006)etc.]. 

We observed that the four dimensional work of M. Sharif and M. Farasat Shamir (2010) regarding non- 
vacuum perfect fluid solutions in the metric version of f (R) gravity can further be extended to the higher five 
dimensional Bianchi type-I space-time and therefore, an attempt has been made in the present paper. 
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The paper is organized as follows : In section-2, we briefly give the five dimensional field equations in 
metric /(«) gravity. Section-3 is used to find exact non-vacuum solutions of Bianchi type-I space-time in v 5 , 
section-4 is dealt with some physical quantities. Section-5 and 6 are dealt with five dimensional models of the 
universe and in the last section-7, we summarize and conclude the results. 

II. FIVE DIMENSIONAL FIELD EQUATIONS IN f(R) THEORY OF GRAVITY 

The five dimensional field equations in / (R ) theory of gravity are given by : 

F(R)R iJ - ~f(R)g v - V, V j . F(R)+g .. DF(R) = kT .. , (i, j = 1,2,3,4,5) (1) 
df (R) 

where F (R ) = , □ = V V . , V . is the covariant derivative and T is the standard matter energy 

dR 

momentum tensor. 

Contracting the above field equations we have 
5 

F(R)R - —f(R) + 4 uF(R) = kT . (2) 
2 

Using this equation in (1), the field equations take the form 

1 

F(R)R -V V F(R)-kT = -[F (R)R □ F (R) - kT ]g . (3) 

In this way we have eliminated f(R) from the field equations and therefore, the equation (3) helps us to solve 
the field equations. 

III. EXACT NON-VACUUM SOLUTIONS OF BIANCHI TYPE - 1 SPACE-TIME IN v 5 

In this section we find exact non vacuum solutions of five dimensional Bianchi type-I space time in f(R) 
theory of gravity. The line element of Bianchi type-I space-time in V 5 is given by 

ds 2 = dt 2 - A 2 (t)dx 2 - B 2 (t)(dy 2 ) - C 2 (t)dz 2 + D 2 (t)du 2 (4) 
where A , B and C are cosmic scale factors. The corresponding Ricci scalar is 

A B C D AB AC AD BC BD CD 

R = -2[—+—+—+ — + + + + + + ], (5) 

A B C D AB AC AD BC BD CD 

where dot means derivative with respect to t . The energy momentum tensor for perfect fluid gives 

Ty = (P + P)UM . - pg .. (6) 

satisfying the equation of state 

p = wp , 0 < ve < 1 , (7) 



where p and p are energy density and pressure of the fluid while u. = J g 00 (1,0,0,0,0) is the five- 
velocity in co-moving co-ordinates. Since the metric (4) depends only on t , equation (3) yields a set of 
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differential equation for F (t), A,S,C,pand p . Thus the subtraction of the 00-component and 11- 
component gives 

B C D AB AC AD AF F k 

-—-—-—+ + + + -— -— (p+ P ) = 0. (8) 

B C D AB AC AD AF F F 

Similarly the subtraction of the 00-component and 22- component gives 

A C D AB B C B D B F F k 

-—-—-—+ + + + -— -— (p + p) = 0. (9) 

A C D AB BC BD BF F F 

The subtraction of the 00-component and 33- component gives 

A B D AC BC CD C F F k 

-—-—-—+ + + + -— -— (p+p) = 0, (10) 

A B D AC BC CD CF F F 

and the subtraction of the 00-component and 44- component gives 

A B C AD B D CD D F F k 

-—---—+ + + + -— -— (p+ p) = 0. (11) 

A B C AD BD CD DF F F 

It is interesting to note that there are three independent differential equations in five dimensional case also. 
The conservation equation, T. ". =0 leads to 

A B C D 

p + (p + p)[ — +—+—+— ] = 0. (12) 
A B C D 

Thus we have four differential equations with six unknowns namely A, B , C , F , p & p . 
From equations (8) - (9), (9) -(10) and (8) - (10), we get respectively 

A B A B CDF 

-—-—+(—-—)+(—+—+—) = 0, (13) 
A B A B CDF 



B C B C A D F 

— -—+(--—)( — +—+— ) = 0, (14) 
B C B C A D F 



CD C D A B F 

— -—+(—-—)(-+—+—) = 0 . (15) 
CD C D A B F 



DA D A B C F 

— -—+(—- — )( — +—+—) = 0 . (16) 
DA D A B C F 



These equations imply that 
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B dt 

— = d,exp[ c.J-j— ], (17) 
A a F 

C dt 

— = d 2 exp[ c 2 f— — ] , (18) 
B a F 

D dt 

— =d 3 exp[ c 3 ;— — ] (19) 
C a F 

A 

— = d 4 exp[ c 4 J^] (20) 
D a F 



where c, , c 2 , c 3 , c 4 and d , , d 2 , d 3 , d 4 are constants of integration which satisfy the relation 

c, + c 2 + c 3 + c 4 = 0 , d { d 2 d 3 d 4 = 1 . (21) 
Using equation (17), (18), (19) and (20), we can write the metric functions explicitly as 

dt 



A 



= ap exp[ q f— — ] , (22) 
J a F 



dt 



B = ap 2 exp[ qt -j— ] , (23) 
J a V 

df 

C = ap } exp[ 4 3 j — — ] , (24) 
a F 

d/ 

D = ap exp[ 9 i — — ] , (25) 
a F 

v , ,-3 ,-2 ,-1.1/4 . , ,-2 .-1.1/4 , , ,2 .-1.1/4 . , .2 ,3.1/4 

where Pj = (d, d, d 3 ) , p 2 = (d l d 2 d 3 ) , p 3 =(d,d 2 d 3 ) , p 4 =(d,d 2 d 3 ) (26) 



and 



3C 1 ~ 2C 2 ~ C 3 C l ~ 2C 2 ~ C 3 £l + 2C 2 ~ C 3 £l + 2C 2 + 3C 3 



c. and d. are constants of integration. Using power law relation between F and a we have 

F =ka m , (28) 
where k is the constant of proportionality, m is any integer ( here taken as -3 ) and a is given by 

a = (nit + k l )"" , n * 0 

a = it 2 exp( It) , n = 0 , (29) 
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where k t and k 2 are constants of integration. It is mentioned here that we have used H = la " , 

I > 0 , n > 0 to get the above equation. Thus we obtain two values of the average scale factor corresponding to 
two different models of the universe. 

IV. SOME IMPORTANT PHYSICAL QUANTITIES 

In this section we define some important physical quantities 

The average scale factor and the volume scale factors are defined respectively as under : 



a = (ABCD ) 4 , V = a = ABCD . (30) 
The generalized mean Hubble parameter H is defined by 
1 

H = -[H l + H 2 + H 3 + //J, (31) 
4 



ABCD 
where H = — , H = — , H = — , H = — are the directional Hubble parameters in the directions 
ABCD 

of x, y, z and u axis respectively. Using equations (30) and (31), we obtain 



1 V 1 a 

H= = -[H l + H 2 + ff 3 + flj = -. (32) 

4 V 4 a 

The mean anisotropy parameter A is given by 

— 1 4 H - H 

A = — z (— ! ) . (33) 

4 < = i H 

The expansion scalar 6 and shear scalar a 2 are defined as under 
A B 2 C 

e = M '=— +— + , (34) 

A B C 



a=-a.a" (35) 

2 " 

1 1 

where ct .. = — (V u . + V u . ) - — g .9 . (36) 
In thermodynamics, the entropy of the universe is given by 

Tds = d(pV)+ pdV . (37) 

V. FIVE DIMENSIONAL MODEL OF THE UNIVERSE WHEN n * 0 

For this model, F becomes F = k(nlt +k l Y il " and the corresponding metric coefficients A, B and C turn 
out to be 
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, , „ qinlt+k')" 

A = pAnlt + k.) exp[ — ! ], n*\ (38) 

kl(n-l) 

rt-l 

(n/f + A:, ) " 

S = p 2 (nfr + kA exp[ — ! ], n * 1 (39) 

kl(n-Y) 

n-i 

. . qAnlt + k') " 

C = p^nlt + k t ) exp[ — ! ], n * 1 (40) 

kl(n - 1) 

,,„ q. (nit + k.) " 

D = p (ntt + *,)"" exp[ ^ ^ ], n*l (41) 

kl(n-l) 

The mean generalized Hubble parameter and the volume scale factor become 

H = , V = (nit + kA 4 '" . (42) 

nit + k 1 

The mean anisotropy parameter A turns out to be 

2 2 2 2 

— q, + + a, + q , 

A = ^ - 1 ^ . (43) 

2k 1 (nit + k x ) ltn 

The expansion 6 and shear scalar a 2 are given by 

a , 2 2 2 2 

41 , q,+q.+q,+ q. 

a 2 =- — — . (44) 



nit + k t 2k' (nit + k.) 2 ' 



For stiff matter (co = 1) , were have p = p . Thus the energy density and pressure of the universe become 



. 2 



- 12 kl q, + q.+q,+q~ 
2kp=2kp = ^ ^ (45) 

/ 7 r 3/n+2 , . , , . 5/n 

(n/f + fej ) & (n/f + ) 

The entropy of universe is given by 

1 3/ 

Tds = -[6kl(2n - 5)(nlt +£ 1 )""~ 3 - (4' + q\ + q\ + q] )( nit + /c 1 )~ 1 "'~ 1 . (46) 

k 2k 

Also equation (12) leads to 

P = ^7, (47) 
V 

where c is an integration constant. It is mentioned here that this value of p , when compared with the value 
obtained in equation (41), gives a constraint 

kc + 6kl 2 = 0 (48) 
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Which holds only when n = 5 / 2 and q 2 + q 2 + q 2 + q 2 = 0 . The function of Ricci scalar, / ( R ) is 
2k 12 

f(R) = (nit + k i y iln R + — kl 2 (2n - 5)(nlt + k l )~ 2 "~ 3 ' n , (49) 

where R = R l = 4/ 2 (2w - 5)( nit + k 1 y 2 . For a special case n = — , f (R) turns out to be 

2 

/(j?) = — (- — -y ,2n -—u 2 (- — -) 3,2n+1 (50) 

5 16/ 5 16 Z 

which gives / (^ ) in terms of R . 

VI. FIVE DIMENSIONAL MODEL OF THE UNIVERSE WHEN n = 0 

Here the metric coefficient take the form 

q , exp( - /f ) 

A = p fc, exp( &)exp[ ], (51) 

klk 2 

a , exp( - // ) 

fi = p 2 k 2 exp( fo)exp[ - ^— ], (52) 

klk , 

q 3 exp( - /f ) 

C = p k exp( fr)exp[ ]. (53) 

klk 2 

a , exp( -It) 

D = p 4 k 2 exp( It ) exp[ - ] . (54) 

klk 2 

The mean generalized Hubble parameter will become 

// = / (55) 
while the volume scale factor turns out to be 

V = k 2 exp( Alt ) . (56) 



The mean anisotropy parameter A becomes 



2 2 2 2 

— q, + q. + q, + q . 

A = t- 2 2 \ ]exp( -2ft) (57) 

4/ fc 2 

while the quantizes 8 and a 2 are given by 

2 2 2 2 

, q, + q.+q,+q. 

9 = 4/, a 2 =[^ ^]exp(-2ft). (58) 

2k kl 



For stiff matter, the energy density and pressure turn out to be 
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2 2 2 2 2 

- 12 kl exp( -3lt ) q, + q. + q, + q , 
2kp = 2kp == -f --[^ ^— ^]exp( -5ft) . (59) 

The corresponding entropy is 

1 31 

Tds = -[-30Jfc* 2 f 3 exp( It)- (q 2 + q\ + q\ + q 2 4 )exp( -It)]. (60) 

k 2 kk 2 

The constraint equation with the condition, q 2 + q \ + q ] + q 2 = 0 , is given by 

kc - 12 kk 5 2 l 2 exp( 5 It ) = 0 . (61) 
The function of Ricci scalar, f(R), takes the form 

2k , 

/(*) = -exp( -3ftX*- 30/ ). (62) 

5k 2 

2 \2l 2 k %n 

which reduces to f(R)=-[ f' 5 (R - 30 1 2 ) (63) 

5 kC 

using the constraint equation (62). This corresponds to the general function / ( R ) 

/(*)= X a n R " . ( 64 ) 
where n may take values from negative or positive. 

VII. CONCLUDING REMARK 

(i) In this paper we have studied the expansion of the universe in metric / (s ) theory of gravity. Using the non- 
vacuum field equations, we have obtained exact solutions of the Bianchi type-I space-time. These exact 
solutions correspond to two models of the universe. For n * 0 , we have obtained singular model of the 
universe and for n = 0 we found a non-singular model. For these solutions we have evaluated some important 
cosmological physical quantities such as expansion scalar 6 , shear scalar a 2 and mean anisotropy parameter 
A . The entropy of the universe is also found. 

Our observations regarding the model of the universe are as under 
(i) For n * 0 , singular model of the universe 

For this model, we have a singularity at t = t = - — . 

nl 

The physical parameters H , , H 2 , H , H 4 H ,9 and a 2 are all infinite at this point for n > 0 but volume 
scale factor vanishes. 

The mean anisotropic parameter A is also infinite at this point for 0 < n < 1 and it will vanishes 
n > 1 . The function of the Ricci scalar f(R), energy density p , pressure p and T are also infinite while 
the metric function A , B and C vanish at this point of singularity. 
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The model suggests that expansion scalar 8 and shear scalar a decrease for n > 0 with the time. 
The mean anisotropic parameter also decreases for n > 1 with the increase in time. 
This indicate that after a large time the expansion will stop completely and universe will achieve isotropy. 

The isotropy condition i.e, a 2 / 6 0 as f -> oo , is also satisfied. The entropy of the universe is 
infinite for n > 1 / 3 . 

Thus we can conclude from these observations that the model starts its expansion from zero volume 
with infinite energy density and pressure at t = t s and it continues to expand with time. 

(ii) For n = 0 , non-singular model of the universe 

The physical parameters H t , H 2 , H , , H 4 , H , a 2 and A are all finite for all finite values of t . 

The mean generalized Hubble parameter H and expansion scalar 9 is constant while / (R ) is also finite here. 

The metric functions A , B and C do not vanish for this model. 

The entropy of the universe is finite. 

The energy density and pressure become infinite as t —> -oo which shows that the universe started its 
evolution in an infinite past with a strong pressure and energy density. 

The isotropy condition is also verified for this model. 

The volume scale factor increases exponentially with time which indicates that the universe starts its expansion 
with zero volume from infinite past. 

We observed that all the four dimensional results of Sharif and Farasat (2010) can be obtained from our 
investigations after the reduction of dimension.. Therefore it is pointed out that the solution of Sharif and Farasat 
(2010) is a particular case of solution presented here. 

We think that this new exact higher dimensional solution should bring some additional information and 
therefore, they need to be further investigated. 
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